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h  new  parametric  SCF  MO  method  (AMI)  has  been  developed  which  seems  to  have 
overcome  the  remaining  deficiencies  of  MINDO/3  and  MHDO,  in  particular  their 
failure  to  reproduce  hydrogen  bonds,  overestimation  of  nonbonded  repulsions, 
and  poor  results  for  compounds  containing  nitrogen  end  oxygen,  specifically 
those  containing  the  N02  group.  MHDO  parameters  have  been  devised  for  bromine 
and  iodine,  d-AOs  have  been  introduced  into  WDO  and  the  resulting  treatment 
(MHDOD)  parameterized  for  P  and  S.  This  work  was  suspended  to  take  advantage 
of  AMI.  The  programs  have  been  rewritten,  improved,  and  combined  into  a 
single  package  (MOPAC)  which  will  be  sent  to  QCPE.  The  new  programs  run 
nearly  three  times  faster. 
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j  ABSTRACT 

X  C:.' 

A  new  parametric  SCF  MO  method  (AMI  has  been  developed  whioh  seems 

to  have  overcome  the  remaining  deficiencies  of  MIND0/3  and  MNDO ,  in 

particular  their  failure  to  reproduce  hydrogen  bonds,  o ve r es t ima t i on  of 

nonbonded  repulsions,  and  poor  results  for  compounds  containing  nitrogen 

and  osygen,  in  particular  those  containing  the  NO  group.  MNDO 
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parameters  have  been  devised  for  bromine  and  iodine.  d-AOs  have  been 
introduced  into  MNDO  and  the  resulting  treatment  (MNDOD)  parameters  for 
P  and  S.  This  work  has  been  suspended  to  take  advantage  of  AMI.  Our 
programs  have  been  rewritten,  improved,  and  combined  into  a  single 

package  (MOPAC)  whioh  we  will  be  sending  to  QCPE .  The  new  programs  run 

nearly  three  times  faster 

INTRODUCTION 

Unespeoted  practioal  problems  Interfered  with  the  project  supported 
by  this  grant  during  the  first  part  of  the  grant  period.  I  had  obtained 
a  grant  from  the  National  Soience  Foundation,  with  matching  funds  from 
the  University  of  Texas,  to  acquire  a  Digital  VAX  11-780  computer  for  my 
theoretical  work.  Due  to  delay  by  the  University  in  providing  their 
contribution,  the  computer  was  not  ordered  for  nearly  a  year  after  I  had 

obtained  the  grant.  Sinoe  the  CDC  6600  computer  at  the  University  of 

Texas  Computuation  Canter  had  become  grossly  overloaded,  my  research 
program  was  brought  to  a  standstill.  Worse  still,  just  before  the  VAX 
was  finally  installed,  the  University  replaced  the  CDC  6600  by  a  CYBER 
750/760  dual  computer.  Since  the  computation  canter,  rather  unwisely, 
wrote  their  own  operating  system  for  this  (as  they  had  done  with  the 
6600)  and  since  it  took  them  several  months  to  get  it  properly 
funotioning,  we  were  left  in  the  situation  of  having  to  rewrite  our 
programs  for  two  different  oomputers  at  the  same  time  without  any 
referenoe  as  an  aid  in  debugging  them.  However,  in  spite  of  these 
obstaoles ,  satisfactory  progress  has  been  made  in  a  number  of  areas. 

DEVELOPMENTS  IN  TECHNIQUE 

The  most  notable  result  during  the  current  grant  period,  as  during 
the  first,  has  been  the  development  of  a  new  parametrio  procedure. 
M1NDO/3  and  MNDO  suffer  from  oertaln  speolfio  weaknesses  that  limit 
their  range  of  applicability,  in  particular  inability  to  reproduee 
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P  a  g  e  3 

hydrogen  ooius  .  ovite;  ii.uitioi  of  nU  [  atonic  :epjIsions  thit  l<;ad  to 

too  positive  energies  for  crowded  molecules  <e  a.  neopentane)  end 

activation  energies  for  reactions  that  tend  to  be  too  large,  and  too 

nooative  energies  for  email  rings.  MNDO  also  gave  poor  results  for 

comoounds  containing  the  NO  croup  which  made  it  unsatisfactory  for 
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studies  of  explosives.  All  those  oroblems  have  been  overcome  in  the  new 
"third  generation"  treatment,  which  we  are  calling  Austin  Molecular 
Model  1  (AMI).  AMI  is  currently  undergoing  final  tests  and  tuning,  to 
make  sure  that  it  is  in  its  optimum  form  before  we  estend  it  to  other 
elements.  The  initial  p a r amo t r i s a t i on  is  for  the  "organic"  elements 
(C ,H,N,0) . 


In  the  meantime  MNDO  has  been  parametrised  for  bromineC63  and 
iodinedl  ,  in  response  to  requests  from  many  organic  chemists,  and  also 
for  aluminumCSl.  Pipers [4 ,51  reporting  results  and  parameters  for 
sulfur  and  chlorine  are  also  in  press.  We  have  also  introduced  d  AOs 
into  MNDO  and  parametrised  the  resulting  treatment  (MNDOD)  for  sulfur 
and  phosphorus.  A  preliminary  pa r ame t r i s a t i on  for  ohromium  was  also 
carried  to  a  point  where  success  seemed  assured.  Further  work  in  this 
area  has  naturally  been  suspended  awaiting  the  finalisation  of  AMI. 

Technical  developments  during  the  current  grant  period  have  led  to 
major  improvements  in  our  computer  programs.  Implementation  of  a  new 
matrix  d i agona 1 i sat i on  a  1  go r i t hml 6 1  ,  due  to  Pulay,  has  shown  it  to  bs 
four  times  faster  than  the  ICEPAK  program.  Together  with 
improvements  in  the  Da v i d on-F 1 e t c he r -Powe 1  I  method  used  for  geometry 
optimisation,  and  various  other  improvements  in  our  computer  programs, 
those  changes  -have  reduced  the  time  required  for  MNDO  calculations  by  a 
faotor  of  thrae.  Wa  hava  also  amalgamatad  our  programs  for  looating 
transition  statss,  for  calculating  foroa  oonstants  and  vibration 

frequanclas,  and  for  estimating  thermodynamic  properties  fr»a  the 
latter,  with  the  basio  MINDO/3  and  MNDO  programs  into  a  single  paokage 
and  in  addition  simplified  the  input.  This  package  should  soon  bs 
available  from  QCPE.  Ve  have  already  supplied  it  to  the  Air  Foroe 

Academy  and  other  DOD  1 aborator tee .  We  have  aleo  Implemented  e  eerlee 
of  graphtos  programs,  written  by  Dr.  David  Pansak  at  du  Pont  de  Nemours 

Co.  and  Capt .  Donald  Storoh  at  ths  Air  Foroa  Aoademy,  for  the  display 

and  plotting  of  molooulas  in  various  formats. 
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APPLICATIONS  TO  CHEMICAL  PROBLEMS 

Th*  r*st  of  th*  work  carried  out  during  th*  ournnt  grant  period 
ha*  b**n  concerned  with  applications  of  our  procedures  to  a  wide  variety 
of  problems  of  soientific  and/or  technological  interest.  A  brief 
summary  follows. 


A.  Initiation  of  Detonation  in  Esplosives 


We  completed  the  study,  initiated  during  the  previous  grant  period, 

of  mechanisms  of  initiation  of  detonation  in  organic  esplosives 

containing  th*  NO  group. 
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(a)  We  confirm*dC71  earlier  work  by  Dr.  J.  Alster,  at  Picatinny 
Arsenal,  whioh  lndioated  that  th*  preferred  mode  of  initiation  in  th* 
detonation  of  methyl  nitrate  <1),  and  so  by  inference  in  other  nitrate 
esplosives,  is  an  unprecedented  insertion  of  an  osygen  atom  from  on* 
molecule  of  1  into  a  CH  bond  of  another;  vis. 


+  // 

H3C-O-N 


ch2-o-no2 

I 

H 


H3CONO  4  HOCH2ONO2 


The  oalculated  activation  parameters,  and  volume  of  lot ivat ion,  indicate 

7 

that  this  reaction  should  be  10  times  faster  than  fission  into 

radicals  < CH  O  +  NO  >  under  the  conditions  of  a  detonation  shook 
3  2 

wav*.  Bimolacular  reactions  had  not  previously  been  considered  as 
possible  candidates  for  initiation  of  detonation.  This  kind  of  reaction 
may  also  be  involved  in  oases  where  nitro  compounds  aot  as  oxidising 
agents . 


(b)  Th*  observed  activation  energies  for  thermal  decomposition  of 
nitromethane  <2>  and  methyl  nitrite  (3)  are  s i gn i f i o ant  1  y  less  than 

m  m 

the  energies  required  to  bring  about  fission  into  radioals  (CH  + 

NO  or  CH  0  ♦  NO).  MINOO/3  e a  1 ou 1  a t i ons C 7 1  indioat*  that  th* 

2  3 

deoompos i t i on  of  2  takes  plaoe  most  easily  via  a  preliminary 

m 


rearrangement  to  3  and  that  th*  easiest  rout*  for  decomposition  of  3 
is  by  intramolecular  elimination  of  HNO;  i.e. 
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A  A 

CH3NO2 — *  CH3ONO ►  CH2O  +  HNO 
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Calculation*  have  also  b**n  can  itd  out  for  th*  nltrocthyUnii  and 
ni t roacet y 1 *n*s .  Re a r r angtmtn t  of  on*  nitro  group  in  *ach  of  th*s*  to 
nitrit*  is  pr*diot*d  to  talc*  plao*  much  sort  easily  than  that  of  2  to 
3.  Th*  rearrangements  are  moreover  facilitated  by  nitro  groups  at  th* 
other  end  of  th*  multiple  bond.  Th*  aotivation  energies  for  th* 
rearrangements  of  t * t r an i t r o* t hy 1 *n*  (4)  and  d i n i t r o ac e t y 1 ene  (5) 

are  consequently  v*ry  small  (<2S  kcal/mol*).  It  therefor*  seems 
unlikely  that  either  will  be  stable  at  room  temperature,  let  alone  at 
higher  temperatures.  This  conclusion  places  obvious  restrictions  on 
possible  synthetio  routes  to  these  potentially  interesting  species. 

This  work  was  reported  at  a  NATO  meeting  at  Corfu,  Creece  in,  July 
1980;  at  a  joint  workshop  sponsored  by  AFOSR ,  ONR ,  and  ARO  at  Berkeley, 
California,  in  January,  1981;  and  at  a  meeting  organised  by  ARO  at 
Hilton  Head,  South  Carolina,  in  April,  1981.  A  fullC71  aocount  is  being 
submitted  for  publication  (JACS). 

B  Properties  of  Nitrocubanes 

Dr.  Alt t*r ,  at  Pioatinny  Arsenal,  has  argued  convincingly  that 
00 t an i t r oouban*  <4)  should  be  th*  "ultimate"  conventional  explosive, 
oombining  high  density  with  a  phenomenal  energy  release  per  gram.  Any 

organio  chemist  would,  however,  have  grave  doubts  concerning  its 

stability. 

Th*  d*oompos i 1 1  on  of  cuban*  (7)  undoubtedly  takes  plaoe  in  a 

stepwise  manner,  formation  of  th*  biradioal  (8)  being  the  rate 

determining  step,  baoaus*  any  prooess  involving  th*  synohronous  breaking 

of  two  or  more  bonds  would  be  "forbidden".  Sinoe  nitro  groups  usually 

seem  to  stabilise  radioals,  introduction  of  nitro  groups  into  7  might 

have  been  expected  to  facilitate  cleavage  to  8  and  hence  the  rate  of 

#* 

its  overall  deoompos i t i on ,  probably  to  the  extent  of  making  it 

dangerously  unstable.  Our  MINDO/3  oaloulations  reproduced  the  expeoted 
meohantsm  for  decomposition  of  7  and  led  to  a  reasonable  estimate  (35 
koal/mols)  of  the  aotivation  energy.  Very  unexpectedly,  however,  they 
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also  Isd  to  t ho  conclusion  that  nitro  groups  should  stabilis*  7. 
Introduction  of  nitro  groups  at  ths  and  of  ths  breaking  CC  bond  raised 
the  motivation  energy  very  signif ioantly,  suggesting  that  6  will 

m 

probably  prove  stable  enough  to  be  of  practical  value.  This  work,  which 
was  also  reported  at  the  meeting  at  Hilton  Head,  has  made  attempts  to 
prepare  6  a  high  priority  project. 

C  Chemical  Storage  of  Energy 

There  is  considerable  interest  in  the  possible  use  of  thermally 
reversible  photochemical  reactions  for  the  storage  of  solar  energy,  the 
pho t oohemi o a  1  reaction  being  endothermic  so  that  its  subequent  reversal 
liberates  the  energy  so  stored.  Host  of  the  work  in  this  area  has 
centered  on  nor  bo r nad i ene  <  9  >  which,  by  triplet  excitation,  is 

•m 

oonverted  to  quad r i o yc 1 ene  (10).  The  reverse  reaction  is  exothermio 
by  >20  kcal/mole.  The  conversion  of  9  to  10  requires,  however, 
light  of  wavelength  <300  na,  which  comprises  only  1%  of  the  solar 
spectrum,  and  only  a  small  part  of  the  inoident  light  energy  is 
recovered  in  the  subsequent  regeneration  of  9.  Attempts  to  improve 
the  performance  of  9  by  introducing  substituents  have  proved  too 

m 

unpromising  to  justify  the  oost  of  further  experimental  work.  The 
methods  we  have  developed  should  prove  especially  valuable  in  this  kind 
of  situation  because  calculations  for  reactions  can  be  carried  out  very 
quickly  and  at  a  minute  fraction  of  the  cost  of  an  experimental  study 
and  because  the  results  given  by  MINDO/3  or  MNDO  are  good  enough  to 
serve  as  a  guide  to  experiment  by  directing  it  into  promising  lines.  An 
extensive  MNDO  study  of  the  effeots  of  substituents  led  to  the 
unexpected  conclusion  that  fluorine  should  exert  an  especially 
beneficial  effeot,  the  tetrafluoro  derivative  (11)  of  9  being 
predioted  to  absorb  at  longer  wavelengths  than  9  while  the 
exothermicity  of  the  reverse  reaction  <  1  2  «+  11)  is  5  0%  greater. 
While  synthesis  of  11  would  present  difficulties,  the  perfluoro 

m  w 

derivative  (12)  of  9  was  predioted  to  be  almost  as  good  and  should 

mm  m 

be  easily  obtainable  from  9  by  addition  of  chlorine,  f 1 uor i na t i on ,  and 

m 

dsoh 1  or i na 1 1  on .  An  aooount  of  this  work  has  been  submitted  to  JACSC83. 


D  Heohanisms  of  Cy o 1 oadd i t i ons 
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MINDO/3  and  MNDO  calculations  havi  been  oarrisd 
of  reverse  chelotropic  reactions  involving  loss  of  CO 
1.  All  cases  but  one  of  the  reactions  were  predicted 
noneoncerted  manner,  breaking 
biradical;  eg. 


0*X 


+  AJ 


t ; 


out  for  a  var  iety 

or  N  see  Table 
2 

to  take  place  in  a 
to  an  intermediate 


Similar  conclusions  were  reached  in  an  earlier  MINDO/3  study  of  the 
Diels-Alder  reaction  between  ethylene  and  butadiene  and  the  same  has  now 
also  been  found  in  a  detailed  investigation,  using  both  MINDO/3  and 
MNDO,  of  the  Diels-Alder  reactions  beteween  maleic  anhydride  (13)  and 
furan  (14a).  2 -me t hy 1 f u r an  (14b),  and  2 , 5 - d lme t hy 1 f ur an  (14c). 

The  calculated  heats  and  entropies  of  activation  for  all  these  reaotions 
were  in  acceptable  agreement  with  eaperiment.  While  there  are  reasons 
for  believing  that  MINDO/3  and  MNDO  exaggerate  the  asymmetry  of 
reactions  of  this  kind  where  two  bonds  are  predicted  to  break  (or  form) 
at  different  stages  in  a  reaction,  the  prediction  that  such  reactions 
are  not  synchronous  seems  very  strong. 


The  only  exception  found  as  yet  is  the  chelotropic  conversion  of 
the  bicyclic  diasene  (IS)  to  bensene  and  nitrogen.  This,  however,  was 
also  exceptional  in  that  the  calculated  activation  energy  was  very  small 
(  ca  2  kc  a  1 /mo  1  e  >  . ' 


Two  papersC9,101  reporting  this  work  have  been  submitted  to  JACS. 


E  A  New  Rule  Restricting  Reaction  Mechanisms. 

The  work  summarized  above  supports  a  conclusion  drawn  from  our 
earlier  very  extensive  studies  of  pericyolio  reactions,  that  reaotions 
of  two-bond  type,  i.e.  reaotions  involving  the  making  or  breaking  of  two 
bonds,  do  not  normally  take  place  in  a  synchronous  manner.  This  new 
rule  oan  be  rationalised  in  terms  of  the  Evans-Pol any  1  treatment  of 
transition  states  and  it  is  consistent  with  all  the  available 
experimental  evidenoe.  The  only  exceptions  are  reactions  involving 
migration  of  hydrogen,  which  are  exceptional  beoause  of  the  exceptional 
propensity  of  hydrogen  to  form  three-center  bonds,  and  reaotions  where 
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tht  aotlvat ion  energy  ot  the  synchronou*  process  is  so  low  that  no 
advantage  can  be  gained  from  staggering  the  bond  at k i ng / br aak ing 
prooesses  (e.g.  the  decomposition  of  IS).  A  full  aceountClil  has  been 
submitted  to  JACS. 

F  Orbital  Isomerism  in  Radicals. 

Some  years  ago  Dewar  and  Na r a y ana swami E 8 1  pointed  out  that  the 

cyclio  carbene  (16),  formed  by  loss  of  N  from  1?,  might  be 

2 

expeoted  to  exist  in  two  singlet  forms  of  comparable  stability,  one  the 
true  carbene  (16)  and  the  other  a  biradical  (18)  formed  by  transfer 
of  a  lone  pair  electron  into  an  empty  IT"  MO.  16  and  18  are  distinct 
species,  not  components  of  a  resonance  hybrid.  because  their  wave 
functions  differ  in  symmetry. 

Skell  and  his  collaborators  C 1 2  3  have  presented  evidence  for 
isomerism  of  this  kind  in  amido  (19)  and  aoylosy  (20)  radicals  where 

W  w  w* 

the  unpaired  electron  can  occupy  either  a  e'  or  a  t  MO.  Ve  have  now 

studiedC133  a  number  of  radicals  of  this  kind,  using  MNDO ,  with  results 

in  full  aocord  with  the  experimental  evidence.  Amido  and  imido  radicals 

are  predicted  to  exist  in  the  'if'  form,  the  isomers  lying  ca  IS 

kcal/mole  higher  in  energy,  whereas  the  «r  forms  of  acyloxy  radicals  are 

predicted  to  be  the  most  stable.  The  f  acyloxy  radicals  are  predicted 

to  be  unsyxuaet r i oal  (of  20),  as  intuition  would  suggest.  An  apparent 

«»  » 

recent  claimC143  to  the  contrary,  based  on  a  "state-of-the-art"  ab 
Inif  to  oalculation,  proves  on  inspection  to  be  meaningless  beoause  the 
radicals  were  assumed  to  be  symmetrical.  Our  calculations  also 
provided  an  intriguing  explanation  of  e.s.r.  evidenocCUl  that  the 
bensoyloxy  radioal  (21)  is  a  radical  but  has  a  symmetrical 

structure.  MINDO/3  predicts  the  and  isomers  of  21  to  have  similar 
energies,  so  i n t er oonver s i on  of  the  two  isomers  via  the 

symmetrical  isomer  could  be  rapid  even  at  the  low  temperatures  used  in 
the  ESR  study. 

C  Photolysis  and  Thermolysis  of  Pe r pr op  1  ono 1 ao t one 
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While  thermolysis  of  derivatives  of  pe  r  p '  op  i  ono  1  a  c  t  one  (22)  gives 

esolusively  carbonyl  compounds  (22  — *  23),  photolysis  also  gives 

osiranes  (of.  24).  These  results  were  i n t e r p r e t ed C 1 3 3  in  terms  of 

orbital  isomerism  (see  F  above)  in  the  intermediate  biradicals 

(of.  25),  it  being  assumed  that  the  thermal  reaction  leads 

to  '7r  r  ad  i  oa  1  s  where  deoar  boiy  1  a  t  i  on  is  (genuinely)  symmetry  forbidden 

whereas  photolysis  leads  to  eicited  ^  radicals  that  lose  CO  easily, 

2 

the  resulting  bir&dicals  (24)  oyelicing  immediately  to  osiranes. 

The  if  radicals  oan,  however,  undergo  a  concerted  rearrangement  and  loss 

of  CO  ,  conversion  of  the  biradical  into  normal  closed  shell  molecules 
2 

acting  as  the  driving  force  for  d e c a r b os y 1  a t i on . 

This  mechanism  is  inconsistent  with  the  esperiment  evidenoe  and 
theoretical  calculations  (see  F>  which  indicate  that  the  ground  states 
of  acylosy  radicals  are  thee  isomers,  the  isomers  being  the  high 
energy  forms.  We  have  now  studied  the  thermolysis  and  photolysis  of  the 
parent  compound,  22,  in  detail,  using  MINDO/3  and  UNDO.  Both 

reactions  start  with  the  formation  of  the  biradical  (25)  as  its 
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gauche  isomer  (27).  This  oan  undergo  either  internal  rotation  to 

the  more  stable  trans  isomer  (2B),  or  loss  of  CO  concerted  with 

-  --  2 

hydroaen  migration  to  form  (23),  or  loss  of  CO  to  form  24  and 

2 

hence  24.  24  can  also  be  formed  in  one  step  from  28  by  concerted 

e 1 imina t i on-oyc 1 i sa t i on .  The  easiest  mode  of  decomposition  of  25  is 

the  second,  accounting  for  the  formation  of  carbonyl  compounds  in  the 

thermal  reaction.  The  formation  of  24  in  the  photolysis  can  be 

attributed  either  to  a  "hot  molecule"  process,  the  estra  energy  allowing 

direot  loss  of  CO  ,  or  alternatively  a  triplet  state  of  25  may  be 

2 

involved,  in  which  case  the  seoond  mode  of  decomposition  would  be 
inhibited  since  it  would  lead  to  triplet  escited  24.  Rearrangement  of 
triplet  27  to  the  more  stable  trans  form  (28)  would  then  lead,  via 
singlet  27,  to  24. 

This  work  has  extensive  implications  oonoerning  pho t ochemi o a  1 
reactions  in  general  and  also  Indicated  the  potential  usefulness  of  our 
prooedures  for  the  study  of  esoited  state  processes.  A  full 

aooountC14]  has  been  submitted  for  publication  (JACS). 

K  Rearrangent  of  Cy o 1 ohes a - 2 , 5 -d i enones . 


The  photochemical  rearrangement*  of  derivative*  of 

cy o l oheia-2 , 5 -d i enone  (29)  to  b  i  oy o 1 o C 3  ,  i  , 0 Ihexenones  (30)  has  been 
studied  in  detail  by  Zimmerman  who  interpreted  them  in  terms  of  an 
initial  cyclixation  of  the  n  -ft  *  e ::  ~  f  t  e  d  state  of  the  dienone  to  a 
biradical  (31)  followed  by  depromotion  of  a  electron  into  the  singly 
occupied  oxygen  AO.  The  resulting  switterion  rearranges  to  the  product 
(30).  Detailed  MINDO/3  calculation*  now  indicate  that  aTT-TT*  excited 
state  is  in  fact  involved,  this  cycliting  to  an  orbital  isomer  of  31 
in  which  the  oxygen  lone  pairs  are  intact .  This  work,  which  has  wide 
implications  in  organic  photochemistry,  is  in  course  of  publication. 

I  The  Isomeric  Benzynes 

Some  years  ago  MINDO/3  calculations  here!  173  led  to  the  unexpected 

conclusions  that  m-bensyne  (32)  should  be  as  stable  as  the  more  usual 

o-bamym  (33),  existing  in  the  cyclized  form  (34),  and  that  a 

similar  bicylio  isomer  (35)  of  p-benxyne  should  be  a  local  minimum  on 

the  C  H  potential  surface.  These  predictions  were  shortly 
64 

afterwards  confirmed  by  ex p e r imen t C 1 8 ,  1 9 3  .  MNDO  c a  1 c u 1  a t i ons C 2 0 1  have 
now  supported  and  extended  these  conclusions  and  also  confirmed  the 
suggest ionC21]  that  rearrangements  of  h ex a - 3 - ene -  1  ,  5 - d i y nes 

(e  g.  36 — >37)  take  place  via  p-benxyne  as  an  intermediate. 

J  Claimed  Rearrangements  of  Phenyl  Cation 

Hehre  et  alCZ21,  and  Spe r an* a C 2 3 3  ,  have  interpreted  certain 
reactions  in  the  gas  phase  in  terms  of  1,2  hydrogen  shifts  in  phenyl 
oation  (38).  MINDO/3  calculations,  however,  predict  this 

rearrangement  to  have  a  very  large  aotivation  energy  (44  koal/mole). 
The  matter  is  of  some  importance  because  analogous  rearrangements  of 
aryl  cations  have  been  postulated  from  time  to  time  in  other  connections 
and  because  MINDO/3  has  so  far  proven  generally  reliable  in  calculations 
for  a  wide  variety  of  o a r b oo a t i ons .  We  have  now  oarried  out  a  detailed 
MINDO/3  studyC273  of  reaotions  which  might  provide  alternative 
explanations  of  the  experimental  data,  finding  them  to  be  viable.  There 
is  therefore  no  reason  to  doubt  the  prediction  that  hydrogen  shifts  in 
aryl  oat  ions  are  in  fact  very  diffioult. 

J  The  C  H  Potential  Surfaoe 
- - -g- e  - - - - 


A  detailed  MINDO/3  ctudyC281  of  the  C  H  potential  surface  Isd 
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to  t he  location  of  no  less  than  nineteen  minima  corresponding  to 
stable  speoies  with  relatively  low  energies,  together  with  the 

transition  states  for  their  i n t e r c on v e r s  i  ons .  No  other  potential 
surface  of  this  degree  of  complesity  has  ever  been  studied  in  such 

detail.  Apart  from  being  something  of  a  tour  de  force  in  the 

theoretical  area,  this  work  is  significant  in  showing  how  rich  and 
comples  potential  surfaces  may  be,  even  when  only  a  few  atoms  are 
involved  Our  results  emphasise  the  uselessness  of  current  ab 

initio  methods  for  studies  of  this  kind,  required  for  the  reliable 
prediction  of  reactions  and  their  mechanisms.  Our  MINDO/3  study  took  oa 
1000  hours  of  central  processor  time  on  a  CDC  6600  computer  Any  ab 
initio  method  of  comparable  performance  would  have  required  ca  1000 
t imes  as  1 ong . 

The  results  are  also  of  general  interest  because  many  of  the  minima 
corresponded  to  molecules  of  very  unusual  types,  implying  the  presence 
of  unusual  kinds  of  mutioenter  bonding.  This  work  should  stimulate 
synthetic  studies  which  might  open  up  new  areas  of  organic  chemistry. 
Our  calculations  also  indicated  that  tetrahedrane  (39)  should  be 
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stable  with  respect  to  d i s s s o c i a t i on  into  acetylene  but  that  it  should 
rearrange  easily  to  o y c 1 obu t a d i ene  (40).  These  conclusions  agree  with 

mr 

the  results  of  recent  ab  initioC262  and  MNDOC273  calculations  and  by 
the  recent  p r epa r a t i on  I  2 8 1  of  the  t e t r a - t -bu t y 1  derivative  of  39  and 
the  observat ion(281  that  it  rearranges  easily  to  the  corresponding 
derivative  of  40. 

Another  conclusion  from  our  calculations  is  that  acetylene  should 
dimerise  easily  to  a  oarbene  (41)  which  could  clearly  act  as  an 
esoellent  initiator  for  polymerisation.  This  provides  a  satisfying 
esplanation  of  the  fact  that  acetylene  detonates  easily,  but  only  under 
pressure  or  as  a  liquid.  Initiation  of  polymerisation  involves  the 
reaction  of  41  with  another  molecule  of  acetylene  and  therefore 
oompetes  with  the  unimolecular  conversion  of  41  to  40  at  higher 

m  w  <v 

oonoent rat  ions  of  the  aoetylene. 
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The  fact  that  t  hi  conventional  strain  energies  of  oyolopropane 
<42>  and  oyolobutan*  (43)  are  almost  the  same,  in  spite  of  the 

w  m  •  % 

enormously  greater  angle  strain  in  the  former,  has  defied  desoription  in 

terms  of  current  organic  ohemical  theory.  It  has  now  been  pointed 

outCSll  that  this  is  an  espeoted  consequence  of  a  fact  long  known  but 

ignored,  i.e.  that  the  resonanoe  integral  between  two  hybrid  AOs  of  a 

given  atom  does  not  vanish,  even  though  the  AOs  are  orthogonal.  The 

CC  «r  bonds  in  a  paraffin  consequently  interact  with  their  neighbors  fust 

like  the  CC  Tf  bonds  in  a  conjugated  polyene.  This  analogy  suggests  that 

such  s i gmacon j uga t i on  in  c y c 1 o p a r a f f i ns  should  have  oonsequenoes 

similar  to  that  of  normal  < p i ) c on j u g a t i on  in  cyclic  polyenes, 

cy c  1  o pa r a f f i ns  thus  eshibiting  aromaticity  if  there  are  ( 4n  + 2 )  electrons 

in  the  CC  bonds  in  the  ring  and  antiaromaticity  if  there  are  4n.  The 

recognition  that  42  is  in  fact  s i gma a r oma t i c ,  being  isoconjugate 

with  ben s ene , ao c oun t s  for  its  unexpected  stability  compared  with  43 

and  also  for  a  number  of  other  anomalies  in  its  behaviour,  for  example 

o 

the  shortness  (1.31  A  )  of  the  CC  bonds  in  it  and  the  appearance  of 
its  proton  NMR  signal  upfield  by  11  ppm  relative  to  paraffinic 
methylene.  S i gmac on j ug a t i on  also  accounts  for  a  number  of  other 
phenomena  that  have  been  difficult  to  explain  in  terms  of  current 
theory,  e.g.  the  the  pyramidal  structure  of  t-butyl  radical  and  the  bent 
structure  of  triplet  carbene.  Other  phenomena  explained  nicely  in  this 
way  include  the  gauche  effeot,  the  anomeric  effect,  and  through-bond 
ooupling  between  orbitals.  A  treatment  of  chelotropic  reactions  in 
terms  of  s l gma o on j uga t i on  also  demonstrates  the  complete  analogy  between 
them  and  other  eye  1 oadd i t i one ,  making  it  unnecessary  to  treat 
chelotropio  reactions  as  an  independent  class.  A  preliminary  aooount  of 
this  work  has  appeared  in  printC29}  while  a  more  detailed  one  is  being 
submitted  to  JACS. 

L  Mutual  Conjugation  in  Radicals  and  Cations 

Thirty  years  ago,  in  the  papersC301  introducing  PMO  theory,  one  of 
the  topios  treated  was  the  possibility  of  interactions  between 
substituents  attached  to  a  common  conjugated  system,  leading  to 
synergism  of  their  substituent  effects.  It  was  shown  that  suoh 
synergism  should  ooour  only  when  one  substituent  is  of  - E  (electron 
donating  • 1 eo t r oma r i c >  and  the  other  of  +E  (eleotron  aooepting 
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eliotiOBeric)  typ*  and  whan  the  substituents  an  attached  to  positions 
appropriately  related  to  one  another.  In  particular,  substituents  in  an 
even  alternant  conjugated  system  must  be  attaohed  to  positions  of 
opposite  parltyCSOl,  one  starred  and  one  unstarred,  while  in  an  odd 
alternant  radical  (the  only  odd  systea  treated)  both  substituents  aust 
be  attached  to  positions  that  are  activeC301.  D i pheny 1 p i or y 1 hy d r a s y 1 
(44)  was  cited  as  an  ezaaple.  The  tera  mutual  conjugation  was 
suggested  for  the  phenomenon  as  a  whole.  Mutual  conjugation  in  even 
systems  had  been  recognised  and  interpreted  in  terms  of  resonance  theory 
but  no  convincing  explanation  can  be  given  in  such  terms  of  autual 
conjugation  in  odd  alternant  systems. 

In  recent  years  the  latter  effect  has  been  "rediscovered"  by  two 

groupsC321  who  have  suggested  alternative  names 

32a  3  2bc 

(me r os t b i 1 i t a t i on ,  captodative  effect  )  for  the 

phenomenon,  giving  an  inevitably  inadequate  resonance  interpretation.  A 
new  and  better  derivation  has  now  been  foundC331  in  terms  of  PMO 
theoryfSOl,  which  moreover  leads  to  the  unexpected  conclusion  that  the 
•ffeot  should  be  even  more  prominent  in  odd  alternant  cations  than  in 
the  corresponding  radicals.  Detailed  MNDO  o a  I cu  I  a t  i  ons  [  3 3  3  have 

confirmed  this  prediction  which  should  have  wide  implications  in  the 
chemistry  of  suoh  species. 

M  Biomimatic  Cy o 1 i x a t i ons . 

Sterols  and  triterpenes  are  now  known  to  be  formed  in  vivo  by 
multiple  cyollsations  of  open  chain  polyenes,  each  ring  being  formed  by 
intramolecular  addition  of  a  carbocation  to  a  C«C  bond,  the  result  being 
another  carbocation  which  in  turn  oan  attack  another  C»C  bond.  This 
kind  of  reaction  has  found  extensive  use  in  synthesis  and  is  now 
oommonly  described  as  biomimetic  cyclisation.  Doubt  remains, 
however,  concerning  its  mechanism.  In  particular,  in  the  case  of 
multiple  oy c 1 i s a t i ons ,  there  has  been  muoh  disoussion  oonoernlng  the 
timing  of  the  various  steps.  The  faot  that  addition  to  eaoh  C - C  bond 
takes  place  s t e r eos pec  i  f  i  o a  1  1  y  trans  has  been  taken  as  evidence  for  a 
fully  synchronous  meohanism.  The  stereochemistry  of  suoh  reaotions  oan, 
however,  be  explained  in  terms  of  a  stepwise  meohanism  if  the  addition 
of  eaoh  oation  to  a  C «C  bond  leads  to  a  IT  oomplexC  34  3  rather  than  a 
olassioal  earbonium  ion.  This  posslbltty  has  been  undervalued  beoause 
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of  the  f  * i lur«  cf  many  orginio  ohtnis  t  s ,  inoluding  tho»»  working  in  this 
arsa,  to  socspt  ths  complex  thsory. 

Vs  havs  now  studied  this  problsn  by  carrying  out  dstailsd 
calculations  for  the  double  oyclication  of  cis  (45)  and  trans 
(44)  2-undeca-d ,  1 0-d i eny 1  cation  to  cis  and  trans  me t hy 1 t e t r a  1 i ny 1 
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cations  (47).  The  reactions  were  found  to  take  place  in  steps, 
via  complexes  (e.g.  48)  as  stable  intermediates,  with  c onsequen t t 3 4 3 

s ter eospec  i  f i o  trans  addition  to  the  C*C  bond.  There  now  seems  little 
doubt  that  similar  complex  mechanisms  will  be  found  for  other 

biomimetic  cycl iiat ions .  A  paptr  C  3S ]  reporting  this  work  has  been 

submitted  to  JACS  . 

N  Bensene  Dication. 

Recent  experimental  work  has  suggested  that  derivatives  of  the 

bensene  dioation  exist  as  "none  1  ass i c a  1 "  pyramidal  structures  (49), 

analogous  to  the  i s oe 1 ec t r on i o  pyramidal  earborane,  C  B  H  (SO). 

2  4  8  ““ 

Detailed  MINDO/3  calculations  have  confirmed  that  the  pyramidal 
structures  for  such  ions  are  minima  on  the  relevant  potential  energy 
surfaces  but  have  also  led  to  the  prediction  that  the  isomerio 

"classical"  structures  (51)  should  be  comparable  in  stability.  A 

Mr  mr 

number  of  ions  of  this  kind,  derived  from  substituted  bensenes,  have 
been  studied,  using  MINDO/3,  and  the  activation  energies  for 
i  n t e r o on v e r s i on  of  the  isomers,  50  and  51  ,  determined.  Another 
problem  of  interest  is  the  multiplicity  of  the  ground  states  of  these 
dioations.  Hexachlorobenxene  dication  It  as  been  reported  to  have  a 
triplet  ground  state.  Our  calculations  ara  in  agreement  with  this 

observation  but  suggest  that  bensene  dication  itself  should  be  a 

singlet  This  work,  which  has  suggested  a  number  of  promising 

experimental  studies,  is  in  course  of  publicationl36) 

0  Reactions  of  Singlet  Dioxygen 

Mr  mt  mt  Mr  mt  mt  mt  mt  mr  Mr  vseeeeee 
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The  meohanism  of  reactions  of  singlet  (  A  >  oxygen  (O  )  with 

g  2 

olefins  remains  the  subjeot  of  controversy.  MINDO/3  c a  1 cu 1  a t i ons t 3 7 1 
suggested  somu  years  ago  that  the  first  step  is  the  formation  of  a 
peroxirane  (e.g.  32),  but  Goddard  et  alC381  have  claimed  on  the  basis 

of  C,VB  calculations  that  biradicals  such  as  S3  are  lower  in  energy 
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than  the  isomeric  peroxiranes  S t • phens on 1 3 V 1  hag  shown  that  suoh 

biradicr.  l‘s  cannot  in  fact  be  involved  in  these  reactions  and  all  the 
experimental  evidence  can  indeed  be  explained  convincingly  in  terms  of 
peroxirane  intermediates  Further  support  for  this  has  now  been 

provided  by  a  detailed  MINDO/3  study  of  the  peroxiranes  that  can  be 
derived  from  singlet  oxygen  and  1 imonene  (54).  The  produots  from  this 

«*  m 

reaction  have  been  analysed  in  detailtdOl,  esper iaent a  1 1 y .  The 
calculations  provide  a  wholly  satisfactory  interpretation  of  the 
eiper lecntal  results.  This  work  is  being  submitted  for  publication  in 
JACS  . 

P  Refutation  of  the  Frontier  Orbital  Method. 

The  frontier  orbital  method  has  attracted  a  considerable  following 
in  reoent  years  in  spite  of  its  obvious  and  familiar  theoretical 
deficiencies,  presumably  on  the  grounds  that  whatever  these  may  be,  it 
nevertheless  works  in  praotioe.  Attsntion  is  drawn  in  a  paper,  about  to 
be  submitted  to  JACS,  of  an  exception  extreme  enough  to  throw  doubt  on 
the  whole  FO  approach. 

All  theoretical  treatments  so  far  tried,  now  including  MNDQ , 
agree  in  concluding  that  the  frontier  orbital  density  in 

1  0  ,  ? -bor as ar ophenan t h r ene  <55>  at  the  9-position  is  the  least,  or 

"»  <*r 

almost  the  least,  at  any  position  and  certainly  the  least  at  any  CH 
group.  According  to  frontier  orbital  theory,  this  should  then  be  the 
least  reactivs  position  in  the  whole  molecule  to  e 1 eo t r oph i  1  i c  attack. 
In  ptact  loeMU  ,  however,  it  ie  the  most  reactive  of  all!  The 
frontier  orbital  density  at  the  seoond  most  reactive  position  (6;  see 
S3)  ie  also  lsss  than  at  a  number  of  other  lese  reaotive  positions. 

The  frontier  orbital  method  is  generally  satisfactory  only  in  the 
oase  of  alternant  hydrooarbons  and  then  only  because  of  a  ooinoidental 
correspondence  between  frontier  orbital  ooeficients  and  localisation 
energies  that  was  pointed  out  many  years  ago.  Sinoe  there  are  no 
problems  that  cannot  be  treated  more  simply  and  more  effectively  In 
terms  of  the  PMO  method  than  in  terms  of  the  frontier  orbital  method  and 
einoe  the  former  is  both  more  reliable  and  better  founded  In  theory.  It 
would  seem  high  time  that  the  whole  frontier  orbital  approach  was 
pensioned  off,  If  not  exterminated 
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Q  Car  bor  an** 

A  very  detailed  MNDO  study  o(  th*  carboranes  has  been  r ep or t *d C  3 9  3  , 
covering  all  the  compounds  that  have  been  prepared  or  postulated  as  well 
as  a  number  of  others.  The  results  were  in  good  general  agreement  with 
experiment  and  led  to  apparently  reasonable  interpretations  of  the 
experimental  data  in  several  places  where  there  has  been  controversy. 

R  Why  Are  Large  Conjugated  Aromatio  Rings  Aromatic? 

It  has  been  taken  for  granted  for  many  years  that  Hueokel's  rules 
apply  to  conjugated  rings  of  all  sites  and  for  the  same  reasons, 
i.e.  the  specific  resonance  stabilisation  in  rings  containing 
<4n+2)  electrons.  Calculations  here  over  th*  years,  by  progressively 

more  sophisticated  parametric  methods,  have,  however,  agreed  in 

incorrectly  predicting  bonds  to  alternate  in  length  even  in  cyclic 
polyenes  of  quite  moderate  lite;  for  example  in  Boeke lhe i d* ' s 
C141annulen*  derivative  (56)  and  even  in  asulene  (57).  W*  have 
assumed  in  the  past  that  this  represented  a  weakness  of  ZDO  procedures. 
However  HaddonC  401  has  now  found  the  same  to  be  true  for  ab  initio 
calculations  using  th*  STO-3G  and  4-31C  basis  sets. 

It  occurred  to  us  that  the  equalisation  of  bond  lengths  might  in 

fact  be  due  not  to  the  aromatioity  of  these  compounds  but  to  an  unusual 

kind  of  electron  correlation.  In  a  cyclic  polyene,  the  number  of 
electrons  is  equal  to  the  number  of  bonds.  There  will  therfore  be  a 
tendency  for  electrons  of  opposite  spins  to  segregate  themselves  in 
alternate  bonds  and  thus  to  make  all  the  bonds  in  th*  ring  equivalent. 
This  cannot  happen  in  an  open  chain  polyene  beoause  th*  number  of  bonds 
is  then  one  less  than  th*  number  of  electrons.  This  would  account 
nicely  for  the  fact  that  bond  lengths  in  open  chain  polyenes  alternate 
strongly  whereas  those  In  aromatio  rings  do  not,  even  when  th*  rings  are 
very  large.  This  pioture  is  of  course  similar  to  Lowdin's  alternant 
orbital  model,  which,  however,  considered  alternation  of  electrons  only 
between  atoms,  not  bonds.  We  have  oheeked  this  suggestion  by  carrying 
out  sp in-unr es t r i o t ed  (UHF)  MNDO  c a  1 o u 1  a t i ons t 4 1 J  for  a  number  of  large 
ring  systems,  with  results  in  very  good  agreement  with  experiment. 
Slnoe  the  alternation  effeot  should  hold  equally  for  4n  and  <4n*2) 
rings,  it  may  well  be  that  th*  larger  4n  polyenes  will  prove  to  have 
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equal  bond  length*  sine*  t h*  effects  of  ant itroeat ioi ty  ere  espeoted  to 
decrease  with  ring  sise. 

S  Other  Topics  Studied. 

A  number  of  other  topics  have  been  studied.  the  i nve * t i g a t i on* 
being  at  various  stages. 

(a)  MHDO  calculations  for  be r y 1 1 i u»- i nd eny 1  and  fluorenyl 

conpl «s  as  C  42 ]  were  carried  out  as  models  for  the  iron  tricarbonyl 
♦ 

complexes,  B*H  and  Fe(CO)  being  analogous  in  their  ability  to 

3 

bring  three  ligands.  The  results  agree  with  experimental  data  for  the 
latter  . 

(b)  A  detailed  c omp ar i s on C 4 3 1  of  the  results  of  MINDO/3 
calculations  for  a  variety  of  open  shell  systems,  using  the  generalised 
SCF  coupling  operator  and  " ha  1 f -e 1 ec t r on "  methods,  led,  as  espeoted,  to 
similar  results,  thus  justifying  use  of  the  "ha  1 f -e 1 ac t r on"  method, 
which  is  simpler. 

<c)  A  program  was  written  for  MNDO  calculations  for  linear 
polymers,  using  the  tight  binding  approximation,  and  used  to  oaloulate 
the  geometry,  vibration  frequencies,  and  electrical  properties  of 
polyethylene  and  po 1 y per deu t er i oe t hy I ene .  The  results  agreed  with 
ex per imen  t  C  4  4  I  . 

<d)  Calculations  for  a  wide  variety  of  haloolef ins  led  to  a 
satisfying  explanation  of  observations  at  Allied  Corp.  concerning  their 
retention  on  glo  oo 1 umns C 4 3 1 . 

(*)  MNDO  calculations  for  the  allyl  and  crotyl  anions,  and  for 
their  BeH  derivatives  (studied  as  models  of  their  lithium  derivatives, 
parameters  for  lithium  being  unavailable),  led  to  results  in  agreement 
with  the  available  experimental  evidenoe  and  provided  satisfactory 
explanations  of  itt461. 

(f)  UNDO  oaloulations  have  bean  used  to  assign  ultraviolet 
phot oe 1  so t r on  <UPE)  spectra  of  phenoxy  radioalC471  and  other  aryloxy 
radicals  and  also  to  show  that  the  UFE  spectrum  measured  for  what  was 
assumed  to  be  pheny lohlorooarbene  (58)  in  fact  corresponded  to  a 

m  m 

mixture  of  the  oh  1  or ooy o 1 ohsp t a t r i enaca r benes ,  eg  *  39,  formed  by 

m  m 

rearrangement  of  39. 
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(g)  UNDO  calculation*  war*  alco  used  to  interpret  tht  up*  sp*ctra 
of  various  d*t ivat iv*i  of  t -bu t y 1 ao • t y 1 *n* C 4 8 3  . 

<h)  UNDO  calculations  of  th*  vibrational  spectra  of  a  variety  of 
boron  compounds  agreed  with  experiment  where  experimental  data  were 
available  and  also  led  to  various  pr *d i c t i ons C 4 9 3 . 

(i)  Ca 1 ou I  a t i on* C SO  3  of  the  ionisation  energies  of  p -qu i nod ime t hane 

<40a)  and  it*  2/5-dimethyl  derivative  <40b>  led  to  a  reassignment  of 
*  *  %  *  *  *» 

their  UPE  spectra,  c on t r ad i o t i ng  a  claim  to  have  observed  "shake-up" 
bands  . 

(j)  Calculations  for  the  reactions  of  oarbon  atoms  with  carbonyl 
compounds  led  to  a  reinterpretation  of  the  experimental  results  reported 
by  Shevlin,  which  was  susequently  confirmed  by  Shevlin  as  a  result  of 
experiments  suggested  by  our  c a  1 cu 1  a t i ons C 5  1  3  . 

(k)  MNDO  oa 1 cu 1  a t i ons C  5 2 , 5 3 3  were  used  to  assign  the  UFE  spectra  of 

pheny 1 phos ph i ne  and  a  number  of  its  derivatives,  with  results  in 
disagreement  with  the  previously  accepted  assignment.  The  new 

assignments  were  confirmed  by  experimental  studies  suggested  by  the 
caloulat ions . 

(l)  Calculations  for  a  number  of  reactions  of  very  simple  systems, 

e.g.  th*  rearrangement  of  vinylidene  (H  C»C : ) )  to  acetylene,  agreed 

2 

gratifyingly  well  with  those  from  high  level  ab  initio 

o a  1 ou 1  a t i ons C  5 6 3 .  In  our  opinion  the  main  chemical  value  of  the  latter 

lies  in  their  possible  provision  of  data  needed  to  calibrate  of 
parametric  procedures  suoh  as  ours.  A  paper  reporting  this  work  is  in 
pr eparat ion . 

( m )  In  th*  same  vein,  MINDO/3  calculations  have  been  oarried  out 

+  ■»■ 

for  th*  ethane  dication  (C  H  ).  Th*  results  agree  with  those 

2  4 

from  a  recent  high  level  ab  initio  calculation,  predicting  th*  most 
stable  form  to  be  th*  diborane  analog,  (41),  and  predloting  an 

At  At 

lsomerio  struotur*  <42)  to  be  a  second  minimum  on  the  potennt ial 

At  At 

surfaoe.  However,  MIND0S3  also  predicted  th*  existence  of  another 
stable  isomer,  i.e.  43. 

At  at 
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(n)  Calculations  for  S  hivi  shown  it  to  ooour  in  two  distinot 

7 

isomerio  forms,  on#  with  a  boat  <  6  4 )  gsomstry  and  ons  a  ohair  (65)  . 
Whils  saoh  isomsr  undorgoss  pstudorotat ion  without  activation,  ths 
isomsrs  ar*  ssparatcd  by  an  apprsciablo  snargy  barrisr.  This  work  was 
carrisd  out  in  collaboration  with  Dr.  E.  Vasssrman,  at  ths  Central 
Research  Laboratories  of  du  Font  de  Nemours  Co. ,  and  will  be  published 
jointly  in  a  paper  reporting  his  experimental  studies  in  this  area. 


! 


I 
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